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ABSTRACT We report the three-dimensional structure of 
osteogenic protein 1 (OP-1, also luiown as bone morphogenetic 
protein 7) to 2.8-A resolution. OP-1 is a member of the 
transforming growth factor (TGF-/3) superfamily of pro- 
teins and is able to induce new bone formation in vivo. 
Members of this superfamily share sequence similarity in 
Iheir C-terminal regions and are implicated in embryonic 
development and adult tissue repair. Our crystal structure 
makes possible the structural comparison between two mem- 
bers of the superfamily. IVe find that although there is 
limited sequence identity between OP-1 and TGF-^2, they 
share a common polypeptide fold. These results establish a 
basis for proposing the OP-l/TGF-^2 fold as the primary 
structural motif for the TGF-/3 superfamily as a whole. 
Detailed comparison of the OP-1 and TGF-^2 structures has 
revealed striking differences that provide insights into how 
these growth factors interact with their receptors. 



Osteogenic protein 1 [OP-1, also called bone morphogenetic 
protein 7(BMP-7)] was originally isolated from bone based on 
its ability to induce new bone formation in vivo (1). Preclinical 
studies in nonhuman primate models have demonstrated that 
OP-1 is effective in restoring large segmental bone defects (2). 
Evidence that OP-1 is synthesized in the kidney (3) and is 
present in the circulation suggests that OP-1 has therapeutic 
potential for the treatment and management of osteoporosis 
and related metabolic bone diseases. In addition, tissue local- 
ization and other preclinical studies suggest that OP-1 has a 
role in the repair and regeneration of urogenital (4), neuronal 
(5), and cardiovascular (6) tissues. Knowledge of the three- 
dimensional structure of OP-1 is essential for understanding its 
mode of action and for providing a basis for the development 
of small molecule therapeutics. 

OP-1 and transforming growth factor pi (TGF-P2) are 2 of 
>30 homologous proteins in the TGF-/3 superfamily. Members 
of this superfamily have diverse biological activities and play 
critical roles in the migration, proliferation, and differentiation 
of mesenchymal cells during embryogenesis and in the repair 
and regeneration of tissues during postfetal life (7). These 
proteins are synthesized as large precursor proteins that un- 
dergo proteolytic processing at RXXR sites to yield mature 
active dimers of disulfide-linked monomers. Each monomer 
contains a conserved C-terminal 7-Cys domain with 20-92% 
sequence identity among superfamily members (7). 

In this paper we report the three-dimensional structure of 
mature OP-1 to 2.8-A resolution.^ The overall OP-1 monomer 
fold is a Greek key motif with approximate dimensions 60 A 
X 20 A X IS A. In large part, the analogy to a left hand for 
TGF-/32 is true for OP-1 as well (8) (Figs. 1 and 2). The overall 
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FIG. 1. Schematic drawing of the OP-1 monomer fold. The OP-1 
cystine knot constitutes the core of the monomer and consists of three 
disulfide bonds; two, Cys-67-Cys-136 and Cys-71-Cys-138, form a ring 
through which the third, Cys-38-Cys-104, passes. The four strands of 
antiparallel 0-sheet, which emanate from the knot, form two finger- 
like projections. An a-helix, located on the opposite end of the knot, 
lies perpendicular to the axis of the two fingers thereby forming the 
heel of the hand. The N terminus that corresponds to the thumb of the 
hand is unresoWed in our electron density map. Unlike TGF-02 (8, 9), 
this N-terrainal region is not stabilized by a disulfide bond. The 
0-sheets are displayed as arrows and labeled from fi\ through 08. The 
a-helix is displayed as a tube and labeled al. Shown in solid thin lines 
are the intrasubunit disulfide bonds that make up the cystine knot. 
Starting from Gln-36, the residues invoWed in regular secondary 
suucture are /31 (Lys-39 to His-41), 02 (1Vr-44 to Ser-46), 03 (Glu-60 
to Ala-63), and 04 (Tyr-65 to Glu-70) in finger 1; 05 (Cys-103 to 
Asn-110), 06 (Ile-112 to Asp-118), 07 (Asn-122 to Tyr-128), and 08 
(Val-132 to His-139) in finger 2; and al (Thr-82 to Ile-94) in the heel. 

architectural similarity between OP-1 and TGF-02 (8, 9) has 
allowed us to construct a structure-based sequence alignment 
(Fig. 3) from which we compared these two structures to begin 
to understand the chemical and structural elements involved in 
determining the specificity of these proteins to their receptors. 

METHODS 

Structure Determination. Crystals of human recombinant 
mature OP-1 were grown by mixing equal volumes of purified 
protein (10, 11), at 10 mg/ml, with 8% saturated ammonium 
sulfate in 50 mM sodium acetate (pH 5.0) (12). The crystals 
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Fig. 2. Stereoview of the OP-1 dimer, viewed along the crystallographic twofold axis from the N-terminal side. The single intermolecular 
disulfide bond linking Cys-103 of each monomer is shown in yellow and lies exactly on this twofold axis. The OP-1 monomer has an accessible 
nonpolar surface area of 4394 A^, and the dimer has an accessible nonpolar surface area of 6831 A^, resultmg in a hidden area upon dimerization 
of 979 A^ per monomer. We similarly calculate values of 4448 A^ for the monomer and 6852 A^ for the dimer of the TGF-P2 structure 2TGI (8) 
and 4412 for the monomer and 6779 A^ for the dimer of the TGF-02 structure ITFG (9). In both cases, the hidden area upon dimerization 
is 1022 A^ per monomer. 



have the symmetry of space group P3221 with unit-cell dimen- 
sions of a = 6 = 99.46 A and c = 42.09 A. One crystal was 
used to collect a complete native data set to 2.8-A resolution 
at 4°C. Two heavy atom derivative data sets were collected at 
4'*C, one from a crystal soaked for 7 days in 0.3 mM uranyl 
nitrate and the other from a crystal soaked for 8 h in 0.5 mM 
sodium gold (III) tetrachloride (12) (Table 1). The native and 
derivative data sets were integrated and reduced with the 
R-AXis-iic software suite (13) and scaled together with the 
CCP4 program ansc (14). Inspection of the Marker sections of 
the difference Patterson map revealed a single uranyl site. The 
position of the single gold site was determined by employing 
cross-Fourier techniques using the uranyl position as the 
phasing site. The heavy atomjc,y,z parameters and occupancy 
were refined with the program teneyck (15). By using these 
two derivatives and their anomalous signals, an initial phase set 



was calculated to 4.0- A resolution with a mean figure of merit 
of 0.72. The phases were improved and extended to 3.5-A 
resolution by cycles of solvent flattening (16) and phase 
combination (17) using the CCF4 (14) crystallographic pack- 
age. A completely interpretable 3.5-A resolution electron 
density map permitted the unambiguous tracing of the 
polypeptide chain and identification of the amino acids from 
Gln-36 to His-139 by using the graphic program o (18). The 
model was refined with the program xplor (19) by using all 
reflections between 10-A and 2.8-A resolution for which Fobs 
> 2.0o<Fobs). There were no water molecules included in the 
refinement. The rms deviation from ideality is 0.02 A for bond 
lengths and 3.2"* for bond angles. The current/? factor is 22.8%. 

Structural Alignment of Sequences and Calcnlation of 
Structure-Based Sequence Identity. The sequences of OP-1 
and TGF-02 were aligned in Fig. 3 according to the corre- 
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FIG. 3. Structure-based sequence alignment of the OP-1 and TGF-02 7-Cys domams. Residues from the heel region invoked in interchain 
contacts in the dimers of OP-1 and TGF-^2 are white type on black bad^ounds, while the residues that they contact in the finger 1 region (upper 
row) and finger 2 region (lower row) of the other chain are blade type on shaded badcgrounds. In both OP-1 and TGF-P2, the same residue positions 
are involved in making interchain contacts. 
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Table 1. OP-1 structure determination parameters 



atatistics 


Native 


uoid 


uranyi 


Data collection 








Completeness to 2.8 A, % 


87 


84 


92 




6.7 


10.0 


6.0 


Phasing and refinement 












12.5 


10.0 


Heavy atom sites, no. 




1 


1 


Isomorphous reflections, no. 




3590 


1863 


Anomalous pairs collected, no. 




2852 


3697 


Mean figure of merit 




0.52 


0.69 


Resolution, A 


2.8 






Water molecules, no. 


0 






R factor, % 


22.8 






rms bond length, A 


0.02 






rms bond angle, degrees 


3.2 







■^Rmcrsie = S (|/ - /, whcrc / IS the intensity measurement 

for symmetry-related reflections and {/} is the mean intensity. *Riu> = 
£ A|fpH Fp\/2 hFpy where Fp and Fph are, respectively, the native 
and derivative structure faaor amplitudes. 

spending regions of local structural identity in the OP-1 and 
TGF-/32 structures. Alignment gaps were positioned in loop 
regions, which is where the local conformational homology of 
the a-carbon traces tends to be the lowest. The structure-based 
alignment of OP-1 and TGF-^2 was then used as a template 
for the alignment of the 7-Cys domain sequences of other 
TGF-0 superfamily members. Alignment gaps were positioned 
in regions that are loops in both the OP-1 and TGF-j32 
structures. Percent identity between pairs of sequences was 
calculated as the number of identical aligned sequence posi- 
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Fio. 4. Protein model of the cystine knot superimposed on the 
corresponding 2Fo - Fc electron density map. The peptide backbone 
for amino acid residues Cys-67 to Cys-71 runs vertically up the right of 
the figure, while the peptide backbone for residues Cys-136 to Cys-138 
runs up the left side. The core disulfide, Cys-38-Cys>104, passes 
through the ring in the middle of the figure wiUi Cys-104 on the near 
left side and Cys-38 on the far right side. Sulfur atoms are depicted in 
green with disulfide bonds absent. The central position on the long side 
of the knot is a Gly residue (Gly-69) that is conserved among cystine 
knot structures (8, 9, 25-27) and the TGF-J3 superfamily sequences. 
Gly is the only amino acid able to adopt the atypical conformation 
necessary to accommodate the packing of this structure (9). This figure 
was rendered by using the program o (18). 
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Fio. 5. {4) Ribbon rendition of the OP-1 dimer displaying the two 
monomer chains, one in green and the other in gold. Residues in the 
region of the proposed receptor binding ridge with solvent-exposed 
side chains are rendered as atomic spheres. The tips of fingers 1 and 
2 of one chain and the loop at the C-tcrminal end of the helix of the 
other chain make up this ridge. Residues that are located at positions 
of variable sequence composition are white while the more conserved 
residues are red. (B and C) The solvent-accessible surfaces of the OP-1 
(B) and TGF-/32 (C) dimers are colored based on their electrostatic 
potential. Surface regions having an electrostatic potential of -3 ArT 
or less are red while surface regions of +3 kT or greater are blue. 
Neutral regions are green or gold to correspond to the backbone 
ribbons shown in A, This figure was rendered by using the program 
INSIGHT 11 (Biosym Technologies, San Diego). 

tions, excluding gaps, normalized to the geometric mean of the 
lengths of the sequences, and multiplied by 100. 

Residue Contacts. In Fig. 3 residues were identified as 
contacting if the distance between the centers of at least one 
nonhydrogen atom from each side chain was less than the sum 
of their van der Waals radii plus 1.1 A. 

Angle and Dihedral Angle Between the Axes of the Helices 
in the Dimer. The helical axis was defined as the line equi- 
distant from the a-carbons in the helical region. A sequence of 
four points is needed to define the dihedral angle between the 
axes of the helices in the dimer. The two inner points were 
chosen to He on the helical axes adjacent to the a-carbon of 
residue His-84 in OP-1 or His-58 in TGF-^2, respectwely. The 
two outer points were chosen to lie on their respective helical 
axes, but their location is arbitrary. To measure the angle 
between the helices, the first two points used to define the 
dihedral angle were translated to superimpose the inner points. 
The resulting three points define the angle. 
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Surfiace Area Measurements. The amount of nonpolar 
surface area lost during the formation of the dimer was 
calculated by using access (version 2.1) with a 1.4-A probe 
(20). Nonpolar surface area is defined as the contribution to 
the total accessible surface from carbon and sulfur atoms. The 
surface area measurement algorithm in access slices the 
structure into 0.25- A slabs perpendicular to the z axis. As a 
consequence, the results are sensitive to the orientation of a 
structure relative to the z axis (20). To minimize this effect, we 
evaluated three perpendicular and one intermediate orienta- 
tions of each structure. The results of these calculations were 
combined by accepting, for each nonpolar atom, the largest 
accessible area measured among the four orientations. The 
values for TGF-P2 reported here were calculated by using 
coordinates from entry 2TGI (8) and entry ITFG (9) obtained 
from the January 1994 release of the Protein Data Bank (21) 
at Brookhaven National Laboratory. 

Electrostatic Calculations. The solution electrostatic po- 
tentials surrounding the OP-1 and TGF-/32 (ITFG) (9) dimers 
were calculated by using DELPHI (22, 23) (Biosym Technolo- 
gies, San Diego). The calculations were performed by using a 
solvent dielectric constant of 80, a solvent radius of 1.4 A, an 
ionic strength of 0.145 M, and an ionic radius of 2.0 A. The 
interior of the protein was modeled by using a dielectric 
constant of 2.0. Formal charges were used and distributed as 
follows: atoms ODl and OD2 of Asp were each charged -0.5, 
atoms OEl and 0E2 of Glu were each charged -0.5, atoms 
NDl and NE2 of His were each charged 0,25, atom NZ of Lys 
was charged +1.0, atoms NHl and NH2 oif Arg were each 
charged + 0.5, and atom OXT of the C-terminal carboxyl group 
was charged —1.0. 

RESULTS 

The Monomer Fold. The conserved 7-Cys domain of the 
TGF-j3 superfamily corresponds in OP-1 and TGF-P2 (8, 9) to 
a ring structure termed the cystine knot (24). Similar structures 
have been observed in several other growth factors that are not 
part of the TGF-p superfamily (25-27). The role of the cystine 
knot in these growth factors is to define the trifold topology of 
the monomer, which allows for structural variation to occur in 
the external loops (24). Fig. 4 shows a representative section 
of the OP-1 electron density map in the region of the cystine 
knot. 

The finger 1 region of OP-1 is an antiparallel /3-sheet 
containing a 13-residue SI loop (Phe-47-Glu-60) (28) (Fig. 1). 
The structural alignment of the OP-1 and TGF-/32 sequences 
in Fig. 3 places two gaps in the SI loop. The first gap represents 
a deletion in OP-1 that aligns with Arg-26 in the a2-helix of 
TGF-/32. This deletion results in a tighter non-a-helical turn in 
OP-1 as compared with TGF-P2. The second gap corresponds 
to the insertion of Gln-53 in OP-1, which has the result of 
directing both Gln-53 and Asp-54 side chains into the solvent. 
By comparison, in the corresponding region of TGF-P2, only 
Lys-31 is in contact with the solvent. These differences in the 
conformation of the SI loop also result in the conserved proline 
(Pro-59) adopting a trans conformation in OP-1 rather than 
cis, as in TGF-/32 (9). The conformation of the SI loop orients 
six nonpolar residues so they can contribute to a solvent- 
inaccessible interface with finger 2. Of these six, four are 
aromatic (Phe-47, Trp-55, Tyr-62, and TVr-65), and two are 
aliphatic (Ile-56 and Ile-57). In all, the conformation of the SI 
loop backbone places five polar residues (Arg-48, Asp-49, 
Gln-53, Asp-54, and Glu-60) in contact with so^ent. The net 
surface charge in this region is -2, whereas it is +2 for TGF-/32 
(Fig. 5). 

The only a-helix in the monomer is located between the 
third and fifth Cys residues (Cys-71 and Cys-104) (Fig. 1). This 
helix extends for 3.5 turns from residues Thr-82 to Ile-94, is 
amphipathic, and contains a number of hydrophobic residues 



that in the dimer make contact with residues from finger 1 and 
finger 2 of the other monomer (Fig. 3). Several hydrophilic 
residues (Thr-82, His-84, and Gln-88) form one wall of an 
internal soWent pocket near the twofold axis of the dimer, 
while others (Asn-83, His-92, and Asn-95) are in contact with 
the external solvent. The conformation of the loop leading 
from the C-terminal end of the helix back to the cystine knot 
is similar in OP-1 and TGF-P2. By comparison, the loop 
located at the N-terminal end of the helix is 3 residues longer 
in OP-1, resulting in a different fold than in TGF-^2. In OP-1 
this loop has an N-linked sugar moiety attached (Asn-80) (R. 
Tucker, personal communication); there is no corresponding 
glycosylation site in TGF-j32. Further, this loop is uncharged 
in OP-1, whereas it is negatively charged in TGF-P2. 

Finger 2 is the second antiparallel 0-sheet in OP-1 (Fig. 1). 
The polypeptide chain reverses direction between segments 06 
and 07 through a 3:5 turn (29) beginning at residue Asp-118 
and ending at residue Asn-122. In contrast, TGF-02 has one 
less residue in this loop and adopts a 2:2 turn (29). Residues 
Arg-129 to Val-132, located between segments 07 and 08, form 
a peptide bridge that crosses over the C-terminal end of strand 
05 and produces a 180^ twist in the finger 2 antiparallel 
0-structure. A similar structure is observed in other cystine 
knot growth factors; however, the peptide bridge length varies 
(30). Within the monomer, finger 2 makes intrachain contacts 
with finger 1 by contributing aromatic residues Tyr-116, Phe- 
117, and TVr-128 and aliphatic residues Val-114, Leu-115, 
Val~123, Met-131, and Val-133 to a solvent-inaccessible inter- 
face. OP-1 and TGF-02 differ by three charges in the region 
of the finger 2 turn; OP-1 has two negative charges while 
TGF-02 has one positive charge. In the region between the 
turn and the peptide bridge, OP-1 has a net charge of +3 while 
TGF-02 is neutral (Fig. 5). 

Dimer Fold, The OP-1 dimer is formed by the association of 
the heel region a-helix of one monomer with the finger regions 
of the other monomer and is stabilized by a single interchain 
disulfide bond invoking Cys- 103 (Figs. 1 and 2). Formation of 
this disulfide bond brings the cystine knot regions close 
together with the monomers oriented such that the molecular 
twofold axis passes between and parallel to the rings of the 
knots and perpendicular to the axes of the fingers, as in 
TGF-02 (8, 9). In addition to the similarity of subunit orien- 
tation, formation of the OP-1 dimer hides nearly equal areas 
of nonpolar surface (Fig, 2). 

A major difference between the OP-1 and TGF-02 dimers 
is the relative orientation of the helices in the heel region. The 
angle between the axes of the helices in the heel region of OP-1 
is 43^ which is 10° larger than that measured for TGF-02. The 
measured dihedral angle between the helices is -20° for OP-1, 
which is 14" more negathre than for TGF-02. Despite these 
differences in helical orientation, the same helix and finger 
residue positions are involved in making interchain contacts, as 
evidenced by the shaded residues in Fig. 3. 

DISCUSSION 

OP'l as a Structural Paradigm for the TGF-0 Superfamily. 

The OP-1 structure presented here confirms the proposal that 
the TGF-02 fold is prototypical for the conserved C-terminal 
region of the TGF-0 superfamily (8, 9). We have aligned the 
sequences of the C-terminal 7-Cys domain of a number of 
TGF-0 superfamily proteins by using the structurally based 
sequence alignment of OP-1 and TGF-02 as a template (Fig. 
3). The sequence of the finger 1 SI loop in OP-1 has higher 
homology with other TGF-0 superfamily members than does 
the corresponding sequence in TGF-02. Specifically, all 
TGF-0 superfamily sequences, except the TGF-0 family and 
Miillerian inhibitory substance, are best aligned with TGF-02 
by using an alignment gap at the same finger 1 position as 
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OP-1. This suggests that the conformation of the OP-1 H loop 
is more characteristic of the superfamily as a whole. 

Based on the structure-based sequence alignment described 
above, we have constructed a matrix of the resulting percent 
identity between family member sequences (Fig. 6). In this set 
of proteins, OP-1 and TGF-02 have the lowest sequence 
identity, only 36%, yet they are structurally very similar, 
suggesting that these molecules have evolved from a common 
ancestor. Secondary structure analyses of the sequences of 
OP-1 and TGF-P2 show that the heel region of OP-1 has a 
propensity to form an a-helix, whereas in TGF-/32 it does not. 
Stability of the helical conformation of the TGF-/32 heel 
region, therefore, requires interaction of this helix with the 
fmger regions. As noted in Fig. 3, the specific residue positions 
involved in this interaction are essentially the same for both of 
these structures. Indeed, residues at these positions vary little 
in their physical and chemical properties for all the proteins 
listed in Fig. 6. We conclude that all of the TGF-jS superfamily 
proteins in Fig. 6 share the OP-l/TGF-p2 structural motif. 

Implications for Receptor Binding Spedlidty. The members of 
the TGF-fi superfamily exert their biological effects by inducing 
the association of specific type I and type II serine/threonine 
kinase receptors. Various members have been shown to bind 
specifically to type II receptors, which then complex with type I 
receptors to initiate signal transduction (40). 

It has been suggested that the receptor binding regions of 
TGF-02 include residues that are both solvent-accessible and 
lie at positions of heterogeneous composition in aligned 
superfamily sequences (9). Divergent structural features in 
both OP-1 and TGF-P2 occur primarily in the external loops 



of finger 1 and finger 2, the loops bordering the helix, and the 
residues in the N-terminai domain preceding the first Cys of 
the cystine knot. These regions are also solvent-accessible. In 
both the OP-1 and TGF-fil dimcr structures, the tip of finger 
2 and the H loop of finger 1 from one chain, and the C-terminal 
end of the a-helix in the heel of the other chain, form a 
contiguous ridge approximately 40 A long and 15 A wide (Fig. 
14). We propose that this ridge contains the primary structural 
features that interact with the receptor and that binding 
specificity between different TGF-fi superfamily members 
derives from conformational and electrostatic variations on 
the surface of this ridge. Differences in the conformation of 
the finger 1 n loop, which constitutes the midsection of the 
ridge, and in the turn at the end of finger 2, which forms one 
end of the ridge, were noted. More significantly, there are 
striking differences in the surface charge of this region in OP-1 
relative to TGF-/32 (Fig. 5 B and C). The ends of the finger 
regions are negatively charged in OP-1 whereas they are 
positively charged in TGF-^2. This results in a net charge of -4 
for the receptor binding ridge of OP-1 versus +3 for TGF-P2. 
Conversely, the p-strand located C-terminal to the turn of 
finger 2 (^7, Fig. 1) is positively charged in OP-1, whereas it 
is negatively charged in TGF-P2 (Fig. 5 B and Q. These 
observations suggest that electrostatic charge plays an impor- 
tant role in the specific interactions of OP-1 and TGF-/32 with 
their receptors and thus provide a molecular rationale for the 
functional differences between these two proteins. 

We are grateful to Prof. Gregory Petsko for his enthusiastic support 
of this project, insightful discussions, and providing unlimited use of 
the crystallographic facilities at Brandeis University. We are also 
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Fig. 6. Identity matrix for the TGF-^ superfamily. Included in the matrix are members that are more than 36% identical to OP-1 in their 7-Cys 
domains. Sequences were positioned in the matrix with those more distant from OP-1 at the bottom and arranged so as to cluster those of high 
identity. Boxes enclose families of sequences having 50% or higher identity with a majority of the other members of the family; gray areas identify 
groups of sequences having 75% or higher identity. Recombinant^ expressed OP/BMP family members known to be able to make bone are denoted 
by a + symbol to the left. Superfamily members whose three-dimensional structures have been determined (OP-1 and TGF-02) are highlighted 
black to the left. Some superfamily members were not included in the matrix because their level of sequence identity with either OP-1 or TGF-^2 
was too low to permit meaningful structure-inferred sequence alignment (i.e.» MuUerian inhibitory substance, MIS). The sequences are referenced 
in Kingslcy (7). except as follows: UNIVIN (31), SCREW (32). BMP-9 (33). BMP-IO (34). GDF-S (35) [also called CDMP-1 (36)]. GDF-6 (35), 
GDF-7 (35), CDMP-2 (36), OP-3 (37), inhibin ^C (38). and GDF-10 (39). Several sequences in the matrix have alternate names: OP-l(BMP-7), 
BMP-2(BMP-2a), BMP-4(BMP-2b), BMP-6(Vgrl), OP-2(BMP-8). 60A(Vgr-D). BMP-3 (osteogenin). GDF-6 (BMP-13), and GDF-7 (BMP-12). 
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